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Abstract The particle distribution function that de- 
scribes two interpenetrating plasma streams is re- 
investigated. It is shown how, based on the Maxwell- 
Boltzmann-Jiittner distribution function that has been 
derived almost a century ago, a counterstreaming distri- 
bution function can be derived that uses velocity space. 
Such is necessary for various analytical calculations and 
numerical simulations that are reliant on velocity coor- 
dinates rather than momentum space. The application 
to the electrostatic two-stream instability illustrates the 
differences caused by the use of the relativistic distri- 
bution function. 
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1 Introduction 



instabilities — relativistic 



Plasma physics — both MHD (magnetohydrodynamics) 
and kinetic theory — are based on the knowledge of a 
distribution function that provides statistics about the 
average velocity direction (bulk flow) and the devia- 
tion from that mean value (known as temperature). 
Whereas MHD calculations are based on the use of a 
Maxwellian velocity distribution, such is not the case 
in kinetic theory, where the distribution function can 
be calculate d using, e. g., the V lasov equation. Kinetic 



theory (e. g. , ISchlickeiseij|2003() is mostly used when di 



lute plasmas are considered that do not satisfy the con- 
dition of frequent binary particle collisions so that no 
Maxwellian velocity distribution is established. 

In various and already historical work, the classic 
Maxwellian velocity distribution has been generalized 
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see, e.g. 



Jiittner 1911a: Jiittner 



The result was the M BJ (Maxwell-Bq l tzmann- Jiittner ) 
distribution function ( Jiittner 1911al : Jiittner 1911bl ). 



which is essentially given through 
/ oc exp( 



(1) 



where a is a temperature-related parameter. 

However, there are cases when the use of a distri- 
bution function in momentum space is (i) not suitable 
for the analytical calculations at hand; and/or (ii) not 
implemented in the numerical (simulation) code. A re- 
cent example is the investigation of LIDAR (Light De- 
tection And Ranging) Thomson scattering systems in 
ITER (Internati onal Thermonuclear Experi mental Re- 
actor) plasmas ( Beausang and Pruntv 20081) . The cal- 
culation was based on a formula teuchinson 1987 ) for 
the scattered power per unit solid per unit angular fre- 
quency, and thus the particle distribution function had 
to be expressed in terms of velocity variables rather 
than momentum variables. 

Another example is the PIC ( Particle-In-Cell) sim- 



ulation co de Tristan (originally iBunemann 



also, e. g.. lSakai. Schlickeiser. and Shuklal 120041 ) . which 



1993. see 



is used for the investigation of plasma instabilities in 
the context of astrophysical scenarios such as the gen- 
eration of magnetic fields at shock wave sites. A basic 
example of such instabilities is the genera lized filamen- 
tation (or, originally, Weibel) instability (|Weibellll959 : 
Friedlll959l) . 

In this short Note, particle distribution functions are 
considered that describe two interpenetrating plasma 
streams. Such distributions are widely used in plasma 
astrophysics, because in the right reference frame, all 
outflow motion into ambient media and shock wave 
sites can be described by coun terstreaming flows (e.g., 
Tautz and Schlickeiser 2005al ). Examples are solar, 
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stellar, and galactic winds, and the interaction of rel- 
ativistic jets such as that from GRBs (gamma-ray 
bursts) and AGNs (active galactic nuclei) with the in- 
terstellar medium. Especially in the latter cases, one 
needs a distribution function that accounts for the rel- 
ativistic effec ts, which can modify the i nstability rates 
significantly ( Schaefer-Rolffs and Tautj 2008 ) . 

In Sec. m it will be shown how the relativistic dis- 
tribution function for counterstreams in momentum 



space, which has been derived (jTautz and Schlickeiser 
2005a[ ) from the MBJ distribution above, can be rewrit- 
ten in terms of velocity variables. In Sec. [3J the appli- 
cation to the two-stream instability will illustrate the 
significant differences between the relativistic and the 
non-relativistic counterstreaming distribution. Finally, 
the results are summarized in Sec. 21 



To transform the distribution function from Eq. ([2]) 
to velocity space, use ise made of the fact 
( Frankel. Hines. and Dewar 19791 ) that 



5„,2 j3„ 



(4) 



where p and v are momentum and velocity components, 
respectively, and where 



1 



7 



(5) 



is the relativistic Lorentz factor. Thus, the transformed 
distribution function in velocity space reads 

/T = C7'cxpJ-a/ 



2 Counterstreaming Distribution 

In momentum space, the two components of the rela- 
tivistic counterstreaming distribution function can be 
constructed intuitively from t he Maxwell-Bo l tzmann- 
Jiittner distribution fuii ction (IJiittneil llQllai IJiittnei 
191 Ibt ISchlickeiseill2004 : 



Schaefe r-Rolffs and Schlick eiseij|2005 
Tautz and Schlickeiser 2005^ as 



/T = Cexp 



1 + (P|| T Po 



-Pi 



(2) 



where a — mc^/{kBT) is proportional to the inverse 
temperature. By integrating over the whole velocity 
space, the normalization constant C can be evaluated 



(|Schlickeiseij|2004[) to 



C 



a 



4:TT{mc)^K2{a) ' 



(3) 



where K2 denotes the modifi ed Bessel function of the 
secon d kind of order two (a. g., Abramowitz and Stegunl 
19741 ). Note that Eq. ^ does not allow for anisotropic 
temperatures (i.e., thermal velocities) in the direc- 
tions parallel and perpendicular to the counterstream. 
There are several ways to implem ent such tem perature 
anistropies i n cov a riant (see, e.g., Yoon 1989t 



Schlickeiser 2004; Tautz and Schlickeiseil 



semi-relativistic (e . g.. | Zaheer and Murtaza 



2005aD 
' 20071 : 



and 



Tautz and Shalchi 20081 ) distribution functions: how- 
ever, there is no general agreement as to which form is 
best suited. However, temperature anisotropics become 
increasingl y negligible for significant co unterstreaming 
velocities ([Tautz and Schlickeiseil 120071) . Therefore, it 
is appropriate especially for relativistic distributions, to 
neglect temperature anistropies at least for the present. 



(6) 



where v\\ ,v±,vq S [— c, c] and where the total velocity, 



(7) 



is also limited by the speed of light, c. Note that the cor- 
rect relavistic addition theorem for velocities has been 
used and that the perpendicular velocity is decreased if 
the counterstreaming velocity, vq, is increased. 

For small velocities and for non-relativistic temper- 
atures, i.e., vo 'ti c and uth ^ c, a series expan- 
sion of both the argument of the exponential func- 
tion and the Bessel function yields the well-known 
non-relativistic version of the distribution function 



(jTautz and Schhckeiseij|2005aj) . namely 



fnr = (wthVTr) ^ exp 



(«|| T fo)^ + vj 



(8) 



where Vth — ^/'^ksT /m = ^J'l/a is the non-relativistic 
thermal speed, which is proportional to the inverse 
square root of the temperature parameter a used for 
the relativistic distribution. 

The total distribution function that describes two 
counterstreaming plasma components is then combined 
as 



/-er + (l-e)/^ 



(9) 



where e £ [0,1] denotes the relative intensity of for- 
ward and backward stream. Both the parameters a 
and the counterstream velocity, vq, can be allowed to 
have different values for the forward and the backward 
stream. In that case, however, vo becomes dependent 
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a) a = 200, V =0.1c 




b) 



a = 50, = 0.3c 



-0.2 -0.1 0.1 0.2 



c) 



a = 5, = 0.3c 





Fig. 1 The counterstreaming distribution function as a function of W||/c for different values of the normahzed temperature 
a and the streaming velocity uq. Shown are the relativistic version (solid line) from Eq. ([5} in comparison to the non- 
relativistic version (dashed line) from Eq. ([8]). Note that, for small temperatures ksT <^ mc and small streaming velocities 
vq <^ c, the two distributions agree (panel a), whereas, even for ksT < rru? , no agreement at all can be seen (panels c 
and d). Note also that, in panel d, the non-relativistic distribution is superelevated by a factor of 50. 



on a (or vice verse), because of the overarching con- 
dition o f vanishing zero-order c urrent that has to be 
fulfilled |Tautz and Lerchelboovt ). 

Furthermore, the integral over the velocity space is 
then transformed as 



d\ 



27r 



duii 



(10) 



so that the distribution function / from Eq. ^ is nor- 
malized to unity. In that context, note that the nor- 
malization factor C from Eq. ([3]), that was derived for 
the distribution function in momentum space, Eq. ([2]), 
remains unchanged throughout the transformation to 
velocity space. 



3 Two-Stream Instability 

To illustrate the use of the relativistic counterstream- 
ing distribution function, Eq. ([2]), the electrostatic 



two-stream instability (see, e.g., ISchlickeisen 12004 : 



Tautz and Schlickeiseij|2005br ) is evaluated both for the 
non-relativistic and the relativistic distribution func- 
tions. Parameters are chosen as vq = 0.5 c and a = 10, 
thus one has Wth = 0.45 c. The relative intensity of 



the two plasma streams is set to be equal to 1/2. For 
the longitudinal mode, the dispersion relation is given 
through 



dp|| 



dpj 



P±P\\ 



7 (fc||U|| 



(11) 



where Wp — yj ^irnq^ / m is the plasma frequency. 

In the case of the non-relativistic distribution func- 
tion, the double integral in Eq. (1111) can be expressed in 
terms of the plasma dis persion f unction or Z function of 
Fried and Contel ( 196ll ) (see also Tautz and Schlickeiser 
2005bl) . whereas such is hardly possible for the relativis- 
tic distribution function. Here, numerical integration 
methods are used instead. 

The numerical solution of the dispersion relation is 
shown in Fig. Obviously, the growth rate, which 
is defined as the (positive) imaginary part of the fre- 
quency uj{k), exceeds is significantly increased if the 
relativistic distribution function is used. Furthermore, 
both the wavenumber of the maximum growth rate, 
^max, and the wavenumber marking the end of the un- 
stable range, fcond, are both larger than for the non- 
relativistic distribution. This can be understood by 
keeping in mind that, from Fig. [Jl one knows that the 
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Fig. 2 The growth rate for the longitudinal mode of the 
two-stream instability, which is the solution of Eq. pip . 
The solid line shows the case of the non-relativistic function 
from Eq. ((HJ, whereas the dots show the relativistic version, 
which uses Eq. The solution of the dispersion relation 
for the relativistic distribution function is, numerically, both 
more delicate and time-consuming and has therefore been 
limited to a few sample points, indicated by the dots. 

anisotropy is much more clear-cut for the relativistic 
distribution. 

4 Summary and Conclusion 

It has been known for a long time that, for relativistic 
temperatures (or, to be more accurate, for relativis- 
tic thermal velocities), the classic Maxwellian distribu- 
tion function becomes increasingly inaccurate. Such 
is escpecially the case for distribution functions that 
describe plasmas streaming with high velocities, as is 
the case for two counterpropagating components. Be- 
cause there exist wide application ranges of such dis- 
tributions, it is both appropriate and necessary to take 
care of a precise basic construction. 

In this short Note, it has been shown how a particle 
distribution function that describes two interpenetrat- 
ing plasma streams can be generalized to relativistic 
velocities. Based on the MJB distribution function, an 
expression has been derived that explicitely uses veloc- 
ity coordinates. Such is advantageous both for analyti- 
cal calculations and for numerical simulations that rely 
on velocity coordinates. 

As an example to illustrate the application of the rel- 
ativistic counterstreaming distribution and, at the same 
time, to demonstrate the both qualitative and quanti- 
tative differences that result from the use of a correct 
relativistic distribution function, the electrostatic two- 
stream instability has been chosen. It has been shown 
that, due to the sharply emphasized anisotropy pro- 
file of the relativistic distribution function, the growth 



rate is significantly higher and the unstable wavenum- 
ber range is more extended as was the case for the non- 
relativistic distribution. 

In future work, the distribution function proposed 
here should be applied mainly to numerical simula- 
tions. The differences between non-relativistic, semi- 
relativistic, and fully relativistic initial distribution 
functions have to be worked out. Only then can one as- 
sess the necessity of a, for analytical calculations rather 
unwieldy, distribution function as that proposed here. 
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